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Context. The FermifLAT collaboration recently reported the detection of starburt galaxies in the high energy y-ray domain, as well as 
radio-loud narrow-line Seyfert 1 objects. 

Aims. Motivated by the presence of sources close to the location of composite starburst/Seyfert 2 galaxies in the first year FermifLAT 
catalogue, we aim at studying high energy y-ray emission from such objects, and at disentangling the emission of starburst and Seyfert 
activity. 

Methods. We analysed 1.6 years of FermifLAT data from NGC 1068 and NGC4945, which count among the brightest Seyfert 2 
galaxies. We search for potential variability of the high energy signal, and derive a spectrum of these sources. We also analyse public 
INTEGRAL IBIS/ISGRI data over the last seven years to derive their hard X-ray spectrum. 

Results. We find an excess of high energy y-rays of 8.3 a and 9.2 a for 1FGL J0242. 7+0007 and 1FGL J1305.4-4928, which are 
found to be consistent with the position of the Seyfert 2 galaxies NGC 1068 and NGC 4945, respectively. The energy spectrum 
of the sources can be described by a power law with a photon index of T = 2.31 ± 0.13 and a flux of FiooMcV-iooGeV =(8.60 ± 
2.27) x lO^ergcm^s- 1 for NGC 1068, while for NGC 4945, we obtain a photon index of T = 2.31 ± 0.10 and a flux of 
^ioomcV-ioogcV =(1.58 ± 0.32) x 10~ u ergcirT 2 s~'. For both sources, we detect no significant variability nor any indication of a 
curvature of the spectrum. While the high energy emission of NGC 4945 is consistent with starburst activity, that of NGC 1068 is an 
order of magnitude above expectations, suggesting dominant emission from the active nucleus. We show that a leptonic scenario can 
account for the multi- wavelength spectral energy distribution of NGC 1068. 

Conclusions. High energy y-ray emission is revealed for the first time in a Seyfert 2 galaxy. If this result is confirmed in other objects, 
new perspectives would be opened up into the GeV band, with the discovery of a new class of high energy y-ray emitters. 
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1. Infoduction 

Somewhat unexpectedly, the FermifLAT collaboration reported 
the dis covery of fo ur radio-loud narrow-line Seyfert 1 galaxies 
et al.ll2009allbh . 



(Abdo 



Caproni et al. 
warped disc, like 



suggesting that these objects, previously 
undetected at high energies, could constitute a new class of high 
energy emitters. Extending this idea to Seyfert 2 galaxies, we 
show here that two such active galactic nuclei (AGN) among the 
closest and brightest in the X-ray sky, NGC 1068 and NGC 4945, 
are detected at high energies using FermifLAT data. We investi- 
gate whether this emission is dominated by AGN or starburst 
activity. 

NCiC 1068 is an archetypal Seyfert 2 galaxy, located at z = 
0.003786, i.e. 14.4 Mpc away, and harbours a h idden Seyfert 1 
core. It was used bv lAntonucci & Miller! (1 19851) to propose the 
AGN unification. Given the proximity of this spiral galaxy, its 
extens on is well observable in visible light, and it is the clos- 
est, as well as one of the brightest, Seyfert 2 galaxy. This source 
exhibi s both AGN and starburst activitie s in its central region 
(e.g. II ester et al.lll987l: iJaffe et al.ll2004l) . A dusty torus lies in 
its cenxal part, dominating the soft X-ray emission by reflection 
from the central nucleus. A circumnuclear starburst region is lo- 
cated c.t ~1 kpc from the core, dominating the infrared emission 
of the broadband spectral energv distribution (SED) (Thronson 
et al.TI?89). 
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(2006) studied the morphology of the inner 
y due to the Bardeen -Petterson effect, using 
VLBA observations. Cot ton et a showed that the emis- 

sion at 43 GHz from the so-called S 1 central radio component is 
more likely domi lated by thermal emi ssion from the hot inner 
region of the obscuring torus, although Gallim ore et al.l (2004) 
pointed out that the North-East radio component could be dom- 
inated by synchrotron emi ssion. This is corroborated by the re- 
sults of lHonig et tlj d2008l) . who argued that the radio emission 
from the core could be dominated by synchrotron or free-free 
emission, while near-IR data are dominated by the thermal emis- 
sion from the dusty torus. Thanks to n ear-IR interferometr ic data 
using the MIDI instrument at the VLT. lRaban et al.l (120091) found 
the torus to be composed of two thermal components at 300 K 
and 800 K. 

In the high energy domain, Ch andra observations of the core 
of NGC 1068 jo ke et al.ll2003l) showed that the X-ray emis- 
sion is due to photoio nisatio n in the extended, clumpy narrow 
line region. iMatt k al.1 (120041) showed that the neutral reflector is 
Compton-thick, u sing XMM-Newton observations. 

Motivated by the presence of a FermifLAT source, 
1FGLJ0242.7+0007, in the region of NGC 1068 in th e 11- 
months FermifLm catalogue (lFGL. IXbdo et al.ll20l0al) , with 
no proposed counterpart in radio nor in y-rays, we analyse here 
1.6 years of data from the FermifLAT instrument, in order to 
better constrain the origin and properties of the y-ray emission 
in this region. 
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Fig. 1. Left: TS map of NGC 1068 between 100 MeV and 100 GeV. The green ellipses show the 68% and the 95% position errors 
from the 1FGL catalogue, the cyan and magenta circles show respectively the position error (at 68% and 95% CL) for the full data 
set with all the events accounted for, and for front events only. The white contours are taken from an optical image from the Digital 
Sky Survey, showing the extent of the Seyfert galaxy. The red boxed points denote the position of the two quasars nearby NGC 1068 
(see text in Section|4]for more details). Right: same as left panel, for NGC 4945. For clarity, we only present here the position error 
circle for all events. 



NGC 4945 is also a Seyfert 2 galaxy at z = 0.001908 exhibit- 
ing starburst activity in it s central region dlwasawa et al][l993t 
Moo rwood & Qlivall994h . Its emission extends up to soft y-rays 
dPetrv et al .1120091) . as observed by INTEGRAL/SPl It is one of 
the brightest hard X-ray AGN (see lltoh et al.ll2008l and refer- 
ences therein, and Ricci et al., in prep.). NGC 4 945 was found 
to be a Compton thick AGN, based o n GINGA dlwasawa et al.l 
Il993h and INTEGRAL observations dBeckmann et al.ll2009h . 
This source was already reported as a high-energy y-ray emit- 
ter by the FermifLA T collaboration in the 11 -months catalogue 
(Abd o et al.ll2010ah . although the authors did not conclude 
whether this high energy emission is due to starburst or AGN 
activity. 

We present here a detailed analysis of 1.6 years of 
FermifLAT data of NGC 4945 to compare the results to those 
of NGC 1068. 



2. Fermi/LAT data analysis 

We present in the following our analysis of FermifLAT data on 
the sources 1FGL J0242.7+0007 and 1FGLJ1305 .4-4928 re- 
ported in the 1FGL catalogue. We will show that these sources 
can be associated to the Seyfert 2 objects NGC 1068 and 
NGC 4945, respectively. 

We analysed ~1.6yr of FermifLAT data, spanning from 
August 4, 2008 to March 15, 2010, from a region of interest 
of 10° in radius around NGC 1068, using the publicly available 
Science ToolQ, and we followed the unbinned likelihood analy- 
sis scheme presented in lAtwood et al.l d2009h . We always used 
the so-called "diffuse class" events, which are the events de- 
tected by the FermifLAT with the highest probability to be y-ray 
photons, and the P6V3 instrument response. 

Using the gtlike tool and assuming a power-law shape for 
the source spectrum, the Test Statistic (TS, iMattox et ai1ll996h 

1 http://fermi.gsfc.nasa.gov/ssc/data/analysis/software/ 



of the likelihood analysis is 68.6, corresponding approximately 
to a 8.3 cr source detection in the 100 MeV- 100 GeV range. The 
corresponding TS map is shown in Fig. Q] The best-fit loca- 
tion of the source using the gtfindsrc tool is 072000 = 2 h 42 m 46 s , 
£72000 = 0°2'14" with an error circle radius of ~6'1 (68% con- 
fidence level, CL), and is fully compatible with the position re- 
ported in the 1FGL catalogue. The maximum photon energy de- 
tected from the source is 20.0 GeV, located at 0'48 from the po- 
sition of NGC 1068. 

Given the energy dependence of the point spread function 
(PSF) of FermifLAT, we performed a second analysis, using only 
front events for which the PSF is narrowefl Front events are 
those converte d in the top layers of the tracker of the LAT in- 
strument (see lAtwood et al.ll2009l for more details). In this latter 
analysis, the TS of the source is 42, still sufficient to derive a 
position of the FermifLAT excess. The best-fit position of the 
source is then a/2000 = 2 h 42 m 49 s , <5/ 2 ooo = -0°0'30", with an 
error circle radius of ~3'A (68% CL), which is only 2'.1 away 
from the nominal position of NGC 1068. This latter result on the 
position is only marginally compatible with the position reported 
in the 1FGL catalogue. 

Given the angular distance between the Fermi source 
1FGLJ0242.7+0007 and NGC 1068, and its optical extension 
of ~6'5, we propose that this Fermi source is actually associated 
with the Seyfert 2 galaxy NGC 1 068 . 

All the sources reported in the FermifLAT 1 1 -months cat- 
alogue dAbdo et al.l l2010ah within a radius of 15° around 
NGC 1068 were included in the likelihood analysis, and mod- 
elled with power-law spectra. We first identify the sources 
not contributing significantly to the likelihood and remove 
them from the model. The spectral parameters of the remain- 
ing sources are left freely varying. For NGC 1068, we obtain 
F tooMeV-iooGeV =(8.60 + 2.27) x 10~ 12 ergcirT 2 s _1 and T - 



2 see http://www-glast.slac.stanford.edu/software/IS/glast_lat_ 
performance .htm 
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1FGL J 1305. 4-4928. The likelihood analysis on NGC4945 re- 
sults in a TS of 85.3 in the lOOMeV-lOOGeV band, equivalent 
to a 9.2 cr detection. Assuming a power-law shape on the source 
energy spectrum, a photon index of F = 2.31 + 0.10 and a flux 
of FiooMev-iooGev =(1.58 ±0.32) x 10~ n ergcrrT 2 s _I are found. 
The highest energy photon detected is 20.7 GeV, located at Q'.22 
from NGC4945. As for NGC 1068, the use of a broken power- 
law or a log parabola did not improve the likelihood, and no 
significant variability was found in the data (see Fig. |2), which 
are statistically consistent with a constant for both sources. 

The analysis of the whole data set gives a position of a 72000 = 
13 h 05 m 33 s , 5/2000 = -49°26'44" with an error circle radius 
of 3'.2 (68% CL), only 1'6 away from the nominal position of 
NGC 4945, while the analysis for the front events only results in 
a position of 072000 = 13 h 05 m 34 s , <5, 2 ooo = -49°26'49" with an 
error circle radius of 3'3 (68% CL). The source position is fully 
compatible with the results reported in the 1FGL catalogue. 

It should also be noted that for both sources, the spectral 
parameters found are fully compatible with the ones reported in 
the 1FGL catalogue. 
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Fig. 2. Top: Light curve of NGC 1068 in the lOOMeV-lOOGeV 
energy band. The time bins are 60 days wide, and the arrows rep- 
resent 95% CL upper limits. Bottom: same as top for NGC 4945. 
In both sources, no significant variability is visible. 



2.31 + 0.13. We also tried to fit the data with a broken power- 
law or a log parabola, but this did not improve the likelihood. 

We performed two other analyses, focussing on the 1— 
100 GeV band, for all events and for front events only. This en- 
ergy band benefits from a better PSF compared to lower energy 
events, despite a much smaller count rate due to the soft spec- 
trum of the source. We obtain consistent results to those reported 
above, including comparable position error on the source. 

We also investigated the potential variability, performing a 
likelihood analysis in different time intervals, for a time bin of 
60 days, letting only the normalisation of the source free to vary. 
As shown in Fig. [2] no significant variability was detected. 

We followed the same procedure for the analysis of 
FermifLAT data in the region of NGC 4945, for the source 



3. INTEGRAL data analysis 

For the extraction of the IBIS/ISGRI spectra, we used all the 
public data obtained by INTEGRAL as of May 2010, for a total of 
564 and 865 pointings (ScWs) on NGC 1068 and NGC 4945, re- 
spectively. The typical exposure of each pointing is (1 -3)x 10 3 s, 
and only ScWs with an effective exposure longer than 200 s were 
kept, spanning times between December 30, 2002 (revolution 
26) and April 7, 2009 (revolution 791). The total exposure is of 
622 ks for NGC 1068, and of 927 ks for NGC 4945. In the large 
17-80 keV band, NGC 1068 and NGC 4945 were detected with 
a significance of 14.0 <x and 138.8 cr, respectively. 

IBIS/ISGRI has a large field of view of 29° x 29° with a spa- 
tial resolution of 12 arcmin. Note that because of the nature of 
coded mask imaging the whole sky image taken by the instru- 
ment has to be considered in the analysis , because all sourc es in 
the field of view contribute to the signal (tCaroli et al.ll 1987h . 

The ISGRI data were reduced using the INTEGRAL Offline 
Scientific Analysis softwar^l version 9.0, publicly released by 
the INTEGRAL Science Data Centre (ICourvoisier et al J 12003). 
The analysis of the IBIS/ISGRI data is based on a cross- 
correlation procedure between the recorded image on the detec- 
tor plane and a decoding array derived from the mask pattern. 
We created mosaic images of all pointings in 10 energy bins, 
and extracted the spectra using mosaic_spec. 

We used the latest detector redistribution matrix files (RMF) 
and calculated the ancillary response functions (ARFs) with a 
weighted average of the 9 available ARFs, based on the number 
of ScWs within the validity time of a particular ARE 

We have also checked the SwiftfBAT dBarfhelmv et alj|2005h 
spectra of the two sources. The 18 month s SwiftfBAT spectra 
have been extracted from the BAT archive (Seg reto et al.l l2010) 
served by the HEAVENS source results archive (Walter et al., in 
prep.). The ISGRI and BAT spectra are in good agreement for 
both sources. 

A simple power-law fit to the combined INTEGRALflSGRl 
and SwiftfBAT data of NGC 1068 using Xspec results in 
^ 2 /d.o.f. = 10.21/l 1. The resulting photon index is F = 
2.08*Qji, and the flux density is Fn-sokev 
1 0" 1 1 erg cm- 2 s" 1 (both 90% CL) . 



1.86!^ x 



http://www.isdc.unige.ch/integral/ 
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Fig. 3. Relationship between SN rate, total gas mass and y-ray 
luminosity of NGC 1068, NGC 4945, NGC 253, M 82, the LMC 
and the Milky Way. 

The hard X-ray spectrum for the combined ISGRI and BAT 
data of NGC 4945 is not well reproduced by a simple power- 
law, and a more complex model is required. A power-law with 
exponential cut-off yields ;^- 2 /d.o.f.=89. 6/9, which is not satis- 
factory, even though the data clearly show a curvature. Instead, 
an absorbed power-law yields ,y 2 /d.o.f.=54.77/9, with a photon 
index of T = 2.03 + 0.04, and an intrinsic absorption column 
N H = (6.66 + 0.45) x 10 24 cirT 2 . Even though the^- 2 is not en- 
tirely satisfactory, we will use the latter spectral model to build 
the broadband SED of NGC 4945, because it provides a better 
description of the data. A detailed analysis of the hard X-ray 
spectral shape is beyond the scope of this paper. 

4. Discussion 

It shall be noted that two quasars are present in the field of 
view of NGC 1068. The quasar SDSS J024230.65-000029.6 
(z = 2.51) lies ~2.'4 from NGC 1068, and another quasar, SDSS 
J024250.98-000031.6 (z = 2.18, a.k.a. QSO 0240-0012), is 
~2'6 away from NGC 1068. However, given their high redshift, 
it is unlikely that one of these objects could be at the origin of 
the high energy emission, alth ough this can not b e completely 
excluded. Indeed according to lAbdo et all (1201 Obh . the highest 
redshift source detected by the Fermi/LXT is a flat spectrum ra- 
dio quasar with a redshift of z = 3.10. Given the similarities of 
the high energy spectra of 1FGL J0242.7+0007 and NGC 4945, 
it appears however more likely that 1FGL J0242. 7+0007 is as- 
sociated with NGC 1068. 

An obvious test to rule out the starburst origin of the high- 
energy emission would be to detect significant y-ray variations 
from one of these objects, as the emission arising from starburst 
activity is expected to be steady. Due to the lack of statistics, no 
conclusions can be drawn about the variability of NGC 1068 or 
NGC 4945 from their light curves (see Fig. [2]). 

Another way to disentangle the starburst or AGN origin of 
the high-energy emission from these galaxies is to compare their 
y-ray luminosity with those of the famous starburst galaxies 



Fig. 4. Spectral energy distribution of NGC 1068, including the 
Fermi/LAT spectrum. The black and red points are archival data 
from the NED, the red ones denote data taken from the central 
region of NGC 1068. For clarity, we only show the INTEGRAL 
IBIS/ISGRI data in blue in the hard X-rays. The EIC model for 
the outflow is shown in blue, and the corresponding SSC emis- 
sion is shown in thin red and magenta lines for first and second 
order components, respectively. The thick red line shows the sum 
of the different emission components from the large outflow. The 
EIC component from the accretion disc is shown in green. 
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Fig. 5. Spectral energy distribution of NGC 4945, including the 
ferm//LAT spectrum (black points). For clarity, we only show 
the INTEGRAL IBIS/ISGRI data in blue in the hard X-rays. 
The model for the EIC component from the accretion disc is 
shown in green. We show in red the data of NGC 253 as taken 
from the NED, with the FermilLKY spectrum published in Abdo 
et al. (|2010cf) as well as the H.E.S.S. flux measurement from 



lAcero et aU (l2009h . for comparison. The two objects have clearly 
very similar SEDs. The luminosity axis on the right is given for 
NGC 4945. 



NGC 253 and M 82, which are also detected i n the very high 
energy domain by H .E.S.S. (lAcero et al.ll2009h and VERITAS 
(lAcciari et al.l l2009). respectively. Both have a luminosity in the 
100MeV-5 GeV band of the o rder of wl0 40 ergs _1 , as detected 
with Fermi (lAbdo et al.ll2010ch . Computing the luminosities of 
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Table 1. Comparison of the starburst properties of NGC 1068 with NGC 4945, NGC 253, M 82, the LMC and the Milky Way. 



Source 


d 


r 




w 




^SN 6 


M eas c 




^5GHz 


^y/^SGHz 




(Mpc) 






(10 39 er 


gs- 1 ) 


(yr 


-') 


(10'M Q ) 


(lO^ergs" 1 ) 


(10 38 ergs-') 




NGC 1068 


14.4 


2.31 ± 


0.13 


170 ± 


32 


0.20 + 


0.08 


4.4 


2.7 


16.6 


102 


NGC 4945 


3.6 


2.30 ± 


0.10 


19.5 + 


2.9 


0.1- 


0.5 


4.2 


1.4 


2.2 


89 


NGC 253 


3.9 


1.95 d 


= 0.4 


7.2 + 


4.7 


0.2 + 


0.1 


2.5 + 0.6 


1.0 


2.0 


36 


M82 


3.6 


2.2 ± 


0.2 


13.0 ± 


5.0 


0.2 + 


0.1 


2.5+0.7 


1.2 


3.0 


43 


LMC 


0.049 






0.041 ± 


0.007 


0.005 + 


0.002 


0.67 ± 0.08 


0.016 


0.017 


24 


Milky Way 








3.2 + 


1.6 


0.02 + 


0.01 


6.5 + 2.0 


0.5 


0.36 


89 



Notes. 

(o) Photon indices and luminosities in the y-ray band from lAbdo et all d2010d) for M 82, NGC 253, the LMC and the Milky Way. The y-ray 
lumi nosities are given in the 100MeV-5 GeV band, and the photon indices T are given above 200 Me V for an easier comparison with the results 

from lAbdo etaLH2010ch . 

< w SN rate estimates fromlWilson & Ufvestadl d!982l) : lBrietz et alj dl994l) : lMannucci et alj d2003l) for NGC 1068, from lLenc & Tingavl d2009l) for 

NGC 4945, and from lAbdo et alj d2010d and references therein) f or M 82, NGC 253 , the LMC and the Milky Way. 

(c) Gas mass estimates from lSageetalJdl990l) for NGC 1068, from lWeifi et~aTJ d2008l) for NGC 4945, and from lAbdo etail d2010d and references 
therein) for M 82, NGC 253, the LMC and the Milky Way. 



Table 2. Model parameters for NGC 1068 and NGC 4945. 



Source 


Component" 


s b 


B (G) 


r b (cm) 


TQQ 


rLnuc (ergs" 1 ) 


R (cm) 


K (cm" 3 ) 


n\ 


»2 


Tbreak 


/max 


NGC 1068 


1 


1.2 


10 4 


2.0 x 10 1S 


130-520 


1.5 x 10 42 


2.2 x 10 2U 


12.5 


2.2 


3.3 


10 4 


10 b 


NGC 1068 


2 


1.2 


1 


10 15 


10 4 


3.0 x 10 41 


5 x 10 14 


10 6 


2.0 






300 


NGC 4945 


2 


1.2 


1 


2.6 x 10 15 


10 4 


5.0 x 10 41 


5 x 10 15 


10 6 


1.5 






160 



Notes. 

(fl) The component 1 refers to the model for the large outflow, while component 2 points to the model for the Seyfert emission in hard X-rays from 
the accretion disc. 



NGC 1068 and NGC 4945 in the same energy band, for compar- 
ison, we obtain 1.7 x 10 41 erg s _1 and 2.0 x 10 40 erg s _1 , respec- 
tively. 

Following lAbdo et aD (1201 Oct) , we compare the supernova 
rate /?sn, the total gas mass M gas and the y-ray luminosities of 
NGC 1068 and NGC 4945 to the ones of NGC 253, M82, the 
Large Magellanic Cloud (LMC) and the Milky Way, as well 
as their infrared and radio luminosities (see Table [TJ. These 
objects are the only extragalactic sources which are not AGN 
known to emit high energy y-rays. Models attributing the y-rays 
to cosmic-ray processes, as expected in a starburst galaxy, de- 
pend on the product M ms R SN Adding NGC 1068 and NGC 4945 
to this picture (see Fig. 2 in Abd o et aTll2010ch . and assuming 
that their y-ray emission is entirely due to starburst activity, con- 
firms the general trend, but we only find a correlation coefficient 
of 0.44 for a linear relationship between M gas /?sN and L y (see 
Fig. [3}, corresponding to a null hypothesis probability as high 
as 38%, although the number of considered sources is obviously 
too small to definitely conclude. 

The y-ray luminosity and the SN rate of NGC 4945 are 
fully consistent with those of NGC 253 and M82, hence even 
though this object is a composite starburst/AGN, its high en- 
ergy emission detected using Fermi could be explained only in 
terms of its starburst activity. Concerning NGC 1068, the situa- 
tion is more complex. Its SN rate is comparable to those of M 82 
and NGC 253, but its radio and y-ray luminosities are higher 
by a factor ~10. This would suggest that its high energy y-ray 
emission is more likely dominated by the central AGN activ- 
ity. Indeed, removing NGC 1068 from the sample of sources, the 
correlation coefficient for the relationship between M gas /?sN and 
the y-ray luminosity becomes 0.95, with a null hypothesis prob- 
ability of only ~1%. This also tends to prove the peculiar role of 
NGC 1068 in this study. This is also strengthened by the fact that 
radio maps of NGC 1068 clearly show a structured jet, on parsec- 



and kiloparsec - scales, modelled by the outfl ow from the central 
AGN (see e.g. iGallimore et al.ll2004 120061) . On the contrary, 
the radio morphology of NGC 4945 shows an extended emis- 
sion consistent with the optical morpho logy tracing the edge-on 
galaxy (see e.g. Hones & McAdam 1992), indicating a starburst 
emission. 

Assuming that the high energy emission of NGC 1068 is in- 
deed due to the AGN activity, and taking a scenario similar to 
what was proposed in Lenai n et alj d200 8) where the jet is mis- 
aligned with respect to the line of sight, we consider here the 
possibility that the outflow in NGC 1068 could be a high energy 
emitter. 

A large, mildly-relativistic zone of the wind-like outflow, at 
a few tenth of parsecs from the core, could emit high-energy y- 
rays through external inv erse Compton process (EIC) (see e.g. 
Begelman & Sikora 1987). At such distances, the infrared pho- 
ton energy density is still high enough to ensure a significant 
emission while being not too important to prevent high optical 
opacity from pair production. At about 100 parsecs from the 
core, the magnetic field strength is also expected to be low, of 
the order of 10 4 G according to the estimate for the equiparti- 
tion magnetic field of the radio component A, located at a pro- 
jected distance of ~350pc from the core of NGC 1068 (Pedlar 
et al. Il983h . No significant short-term variability is expected 
from such a large emitting zone. 

The stationary emission of a blob of plasma is modelled here 
through leptonic processes. The macrophysics of the blob is de- 
scribed by its bulk Doppler factor 5h, radius rb, and strength of 
the tangled magnetic field B filling it up. The energy distribu- 
tion of the leptons is described by a broken power-law with the 
density normalisation K, the indices n\ and rt2, and the Lorentz 
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factors of the individual particles y n 
lows: 



ybreak and y max , as fol- 



Kv~"' 

Ne(7) = \J„ 2 - ni . 



if 7min < T < 7break 
if ybreak < J < Tmax 



[cm- 3 ] (1) 



The se leptons emit through synchro tron self-Compton (SSC, 
see e.g. Ginzburg & Svrovatskii 1965) and EIC processes. The 
latter is described by the fraction of th e core luminosity TL mw . re - 
processed by the emitting source (see llnoue & Takaharal fl996). 
the temperature T of the seed radiation field described as a black- 
body, and the distance R between the external radi ation field a nd 
the re-emitting blob. Second order SSC emission (iReesll 19671) is 
also accoun ted for in this model. Mo r e details on the mode l can 
be found in | Katarzvriski et"aTI (l200lh . lLenain et al . (2008), and 
iLenainl d2009l) . We use y^r, = 1 in the following. 

In Fig. [4] we present such a model with synchrotron emission 
responsible for the radio emission, while the Fermi/LAT data are 
interpreted as EIC emission with the infrared emission providing 
the seed photons, from a multi-temperature blackbody. We also 
show the contribution from the SSC process, which is negligible 
compared to the EIC emission at the highest energies. We note 
that the contribution of second order SSC is negligible in the 
case of our interpretation of the SED of NGC 1068. The model 
parameters are summarised in Table [2] 

The hard X-ray spectrum observed with INTEGRAL does 
not seem to fit in this picture and we propose that this emis- 
sion originates from EIC processes on another population of lep- 
tons, within hot plasma located in the vicinity of the accretion 
disc. The seed photons would then originate from the accretion 
disc itself as is usually invoked to explain the X-ray emission 
of Seyfert galaxies. Evidence for an accretion disc comes from 
the fact that the soft X-ray spectrum of NGC 1068, e.g. as ob- 
served with XMM- Newton, is dominated by thermal reflection 
emiss ion (see e.g. iKinkhabwala et al. 2002; Pou nds & Vaugharj 
2006). We checked that the synchrotron and the SSC emission 
from this component is negligible compared to the one from the 
large outflow described above. 

It should be noted that imaging atmospheric Cerenkov tele- 
scopes have the ability to strongly constrain the highest energy 
part of the particle energy distribution. If the maximum energy 
of the leptons is y max > 5 x 10 6 , a significant signal could be 
detected from NGC 1068 with H.E.S.S., VERI TAS, MAGIC, or 
the fu ture Cerenkov Telescope Array (CTA) dCTA consortium! 
I2010h . 

An alternative solution from our model would consist in in- 
terpreting the hard X-ray emission from INTEGRAL and the high 
energy y-ray emission from Fermi/LAT as originating from the 
same spectral component. In this case, the overall high energy 
emission would be due to EIC process from the same lepton 
population. The spectral shapes of both INTEGRAL and Fermi 
data then strongly constrain the parameters on the lepton energy 
distribution, for which we obtain n\ = 2.0, «2 = 3.6, ybreak = 700 
and K = 450crrT 3 . However, such a particle energy distribution 
can not account correctly for the SED observed in the radio do- 
main. 

In Fig. we show the SED of NGC 4945, including the 
Fermi/LAT and the INTEGRAL spectra we derived in Section [2] 
and[3] respectively. For comparison, we also show the data of the 
starburst galaxy NGC 253, as taken from t he NED, along wit h 
the Fermi/LAT spectrum as published in Abdo et al] (l2010ch . 
Clearly, the two objects have very similar broadband SEDs, 
which match almost perfectly, strengthening the idea that the 



high energy y-ray emission from NGC 4945 could be well ex- 
plained only in terms of starburst activity. NGC 4945 should 
then be detectable in the very high energy domain, and espe- 
cially with H.E.S.S. given its coordinates, at a flux level similar 
to NGC 253. However, one difference between NGC 4945 and 
NGC 253 resides in their hard X-ray emission. While NGC 253 
is not thought to harbour a central AGN, NGC 4945 exhibits a 
characteristic Seyfert emission, which is modelled as above for 
NGC 1068. The corresponding parameters are also included in 
Table |2] 

5. Conclusion 

We reported the detection of the Seyfert 2 galaxy NGC 1068 
in the high energy (lOOMeV-lOOGeV) domain with the 
Fermi/LAT, and compared the analysis of this source with the 
one of NGC 4945 obtained by extending the observing period 
to 1.6 years. The object 1FGLJ0242. 7+0007 was reported as a 
y-ray source without association in the 1FGL catalogue. We can 
now associate it quite firmly to NGC 1068 with an overall detec- 
tion significance of 8.4 cr. The data do not show any significant 
variability over the whole period. 

The high energy y-ray spectrum of NGC 1068 is consistent 
with a power-law of index T - 2.31 +0.13 and a flux density 
of (8.60 ± 2.27) x 10~ 12 erg cirT 2 s -1 . Compared to the y-ray lu- 
minosities of M 82 and NGC 253, whose high-energy emission 
is dominated by starburst activity, we find a too high y-ray lu- 
minosity of NGC 1068 to be explained only by starburst activity. 
We thus propose a leptonic scenario to interpret the high energy 
emission in terms of external inverse Compton process from an 
outflowing relativistic wind launched by the central AGN. 

NGC 4945 is detected at a 9.2 cr level in Fermi/LAT data, and 
its y-ray spectrum is best described by a power-law of index T = 
2.31 ±0.10 and a flux density of (1 .58±0.32)xl0" n erg cnr 2 s _1 . 
This object has a very similar multi-wavelength SED compared 
to NGC 253, and based on its radio, infrared and y-ray luminosi- 
ties, its high energy y-ray emission is most likely due to starburst 
activity. 

If high energy y-ray emission due to AGN activity is con- 
firmed in other Seyfert 2 galaxies, this would mark the discov- 
ery of yet a new class of high-energy y-ray emitters. New data 
from Fermi/LAT in the coming years will be extremely valuable 
in this regard. 
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